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Abstract. We study the phase diagram of an imbalanced two-component Fermi 
gas in optical lattices of 1-3 dimensions, considering the possibilities of the FFLO, 
Sarma/breached pair, BCS and normal states as well as phase separation, at finite and 
zero temperatures. In particular, phase diagrams with respect to average chemical 
potential and the chemical potential difference of the two components are considered, 
because this gives the essential information about the shell structures of phases that 
will occur in presence of an additional (harmonic) confinement. These phase diagrams 
in 1, 2 and 3 dimensions show in a striking way the effect of Van Hove singularities 
on the FFLO state. Although we focus on population imbalanced gases, the results 
are relevant also for the (effective) mass imbalanced case. We demonstrate by LDA 
calculations that various shell structures such as normal-FFLO-BCS-FFLO-normal, or 
FFLO-normal, are possible in presence of a background harmonic trap. The phases 
are reflected in noise correlations: especially in ID the unpaired atoms leave a clear 
signature of the FFLO state as a zero-correlation area ("breach") within the Fermi 
sea. This strong signature occurs both for a ID lattice as well as for a ID continuum. 
We also discuss the effect of Hartree energies and the Gorkov correction on the phase 
diagrams. 
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1. Introduction 

A major experimental breakthrough in the study of ultracold Fermi gases was the 
realization of spin-density imbalanced or polarized Fermi gases [D [21 El IH El El E] • These 
experiments are believed to shed light also on the long standing question of the nature 
of high- Tc superconductivity [El [9] . In a normal superconductor, the main mechanism 
for pairing is the BCS paradigm, however it is not applicable routinely in the high- 
Tc systems [lOj. Imbalanced Fermi gases allow the study of states with more exotic 
pairing, such as the Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) [HKT^llIS] phase and the 
Sarma or breached pair (BP) phase [HI [15]. These concepts have also been considered 
in other fields of physics, such as condensed matter, high energy and nuclear physics 

IT7] . It has been shown experimentally that the imbalanced gas will exhibit phase 
separation by forming a core of BCS superfluid inside a shell of gas in the normal state 
in a harmonically trapped system. FFLO-type features have been predicted to occur 
in harmonically trapped Fermi gases [iHl [iHl [20] as an interface effect [211 122] • Optical 
lattices are extremely promising in this context, since the lattice enhances the FFLO- 
type pairing due to nesting of the Fermi surfaces p3]. Moreover, optical lattices allow 
to manipulate the effective dimensionality of the system as well as the mobility of the 
particles compared to the strength at which they interact 

Density-density correlations have been used as an indicator for different phases in 
optical lattices for bosonic atoms [23] and the idea is promising also in the fermionic 
case [26]. A density-density correlation tells how strongly the atomic densities at 
different positions are correlated. The reason why density-density correlations can 
be a useful way to measure different phases in optical lattices is that while densities 
can be very similar for different phases, the density-density correlations can still be 
markedly different. As an example, one can mention the Bose- Einstein condensate 
and the Mott insulator phases, which show dramatically different density-density 
correlations [23 [2S] • In the Mott insulator one can see clear correlation peaks in the 
density-density correlation, but the density-density correlations vanish (after subtracting 
the product of average densities) for a Bose-Einstein condensate. Likewise, it is possible 
to detect pairing effects in the Fermi gas by measuring density-density correlations 
between different components, as well as antibunching in the density-density correlations 
in a single component Fermi gas |28j . 

In this paper, we study phase diagrams for polarized Fermi gases in optical lattices, 
taking into account the following states: BCS, BP, FFLO, normal state, and phase 
separation into normal and BCS regions. We account for the effect of the harmonic 
trap with the use of the local density approximation and show how shell structures 
such as FFLO - normal state will appear in the trap. We also show how the effective 
dimensionality of the lattice (3D, 2D or ID) affects the phase diagrams and explain 
the relation to the Van Hove singularities of the lattices. Furthermore, we study these 
states by measuring the density-density correlation function of the system. We show that 
from the structure of the the correlation peaks one can clearly distinguish between the 
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spatially modulated FFLO-states, the usual BCS-state, and gain valuable information on 
the structure of quasi-particle dispersions. Finally we discuss the effects of the Hartree 
and Gorkov corrections on the phase diagrams. We focus on the density imbalanced case, 
but the results, especially the phase diagrams for fixed chemical potential difference, 
are also relevant for the mass imbalanced case. Different masses could be introduced as 
different effective masses originating from different hopping strengths, or by considering 
mixtures of fermions of non-equal mass [22] ■ 

Note that by 1 and 2 dimensional lattices we mean here: ID gas in a ID lattice, and 
2D gas in a 2D lattice. The first one could be realized by a 3D optical lattice where the 
confinement is very strong in two directions and intermediate in the third one, and the 
second by strong confinement in one direction and intermediate in two. This is different 
from what is often meant by 1 and 2 dimensional optical lattices superimposed on a 3D 
gas, namely that the ID lattice forms 2D "pancakes" and 2D lattice forms ID "tubes". 
In other words, we study here the actual dimensionality of the lattice, not the effective 
dimensionality of homogeneous space produced by a lattice. 

2. Mean field attractive Hubbard model 

We consider the mean field attractive Hubbard model in the lattice, 

(^Tk4kCTk + ^ikclk^ik + Ac{k+q4-k+q + Aq_k+qCTk+qj ^ Jj , (l) 

k 

where we have limited the study to include only the lowest energy eigenstate of each 
lattice site, i.e. the lowest band. Here the single particle dispersion is 

^ak = 2/^(1 - cos/c^.) + 2Jy{l - cos ky) + 2J^{1 - cos/cj - /i<^. (2) 

The interaction term Acj.k+qC-fk+q corresponds to a plane wave (FFLO) ansatz for the 
order parameter: ?7(c|xC|x) = Ae^*'^'''. The interaction and hopping parameters, U and 
Jj, are defined as in [SU]. Throughout the paper, J without an index stands for the 
largest Jj. Diagonalizing the Hamiltonian with the standard Bogoliubov transformation 
yields the quasiparticle energies 

i?±,k,q = ^^"^'^ ; ^^"-'^ ± ,^{^^^^^^^7^^ (3) 

and the grand potential of the system is 

^ = + E {^l-^^- + ^-.M - ^ 1- ((1 + -"'^^'''■'^) (1 + e^^--))) , (4) 
where {3 = 1/kBT. 

We map out the phase diagrams of the system in section [S] by minimizing Q with 
respect to A and q by keeping the average chemical potential fiaver = (/^t + and 
the difference between the chemical potentials, 6fi = fi-^—fii, fixed. We have also studied 
phase diagrams for fixed filling factors and ui , in section [6l In this case the relevant 
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thermodynamical quantity to minimize is the Helmholtz free energy, which is related to 
the grand potential hj F = Q + fi-^n-i + fi[ni. The major difference between these two 
schemes is that with fixed chemical potentials the BCS state cannot support a finite 
polarization at zero temperature, and BP is not stable, but with fixed densities BCS is 
essentially a special case of the BP state and can have different numbers of the different 
spin components. See e.g. [31j for a more detailed discussion on the stability of different 
phases in these two situations. 

3. Phase diagrams 

Changing the difference of the chemical potentials, Sfi = fi^ — fii, changes the stable 
phase. At T = 0, when 6fi = 0, the result is the standard BCS-state, with A > and 
q = and with equal densities in the (pseudo) spin components. As 5/i is increased 
beyond the Clogston limit [32], \/2Ao, where Aq is the energy gap at 5/1 = 0, the gas 
switches to a polarized FFLO state, with some finite q, which is found by minimizing 
(jlj). This transition is of the first order. As 6fi is increased further, A approaches zero 
and |q| grows, until the gas undergoes a second order phase transition to a polarized 
normal gas. 

Some typical phase diagrams as functions of fiaver = (/^t + and Sfi are given 

in figures [TE In some of the diagrams, some of the circles representing the FFLO - 
normal state phase boundary have fallen onto the boundary between FFLO and BCS. 
This is an artefact of the numerical method employed in the calculations. 

3.1. Van Hove singularity and dimensionality 

The phase diagrams [HIS] show interesting behaviour regarding the shape of the phase 
boundary between the FFLO and normal states. The reason why FFLO can sustain a 
finite polarization, i.e. Fermi surfaces of unequal size, and still be favorable compared 
to the normal state, is that it allows the Fermi surfaces to be partially matched. In 
FFLO, Cooper pairs have a finite momentum 2q, which effectively means a relative 
displacement of the Fermi surfaces by an amount q in order to match them. In 3D and 
2D this nesting of the Fermi surfaces is optimal around the Van Hove singularity, i.e. 
when the Fermi surfaces touch the edge of the first Brillouin zone, because the shape of 
the surface at such densities is octahedral in 3D and square in 2D. In the non-interacting 
case this happens for component a when /ig- = 4 J. Displacing two octahedra or squares 
so that their corners connect actually allows to connect an area, optimally half of the 
minority Fermi surface, instead of just one point, as in the case of spherical Fermi 
surfaces. For this reason the area occupied by FFLO in the phase diagram is more 
dominant than observed for spherically symmetric systems. Moreover, as our results 
show, the Van Hove singularities lead to striking features in the FFLO-normal state 
phase boundary. 

The phase boundary between FFLO and normal state shows special features at the 
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Figure 1. The phase diagram of a two-component Fermi gas in a 3D lattice, at zero 
temperatm'e, as function of the average chemical potential fj-aven and the difference 

(5/i. The dashed line ( ) shows where iif — 4J, that is, the place of the Van 

Hove singularity for the majority component. The dash-dotted line ( — • — ) shows 
where /i| = 4J, which is the minority component singularity. Here J is the hopping 
amplitude. The solid line shows the calculated data points for the BCS-FFLO (in some 
points BCS-normal gas) phase boundary and the circles show the FFLO-normal gas 

phase boundary. The dotted ( ) horizontal lines correspond to the shell structures 

in figure [6l The phase diagram was obtained by minimizing the grand potential ([4]) in 
each point of the diagram. 
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Figure 2. The phase diagram of a two-component Fermi gas in an effectively 2D 

lattice, at zero temperature. The dashed line ( ) shows where = 4J, the Van 

Hove singularity for the majority component. The calculated data points at the phase 
boundaries are marked so that solid line is BCS-FFLO and circles represent FFLO- 
normal gas. 




Figure 3. The phase diagram of a Fermi gas in an effectively ID lattice, in zero 

temperature. The dashed line ( ) shows where /i| = 4J, the Van Hove singularity 

for the majority component. The calculated phase boundaries are shown so that solid 
line is BCS-FFLO and circles represent FFLO-normal gas. 



points where the chemical potential of the majority component, /i| = fiaver + Sfi/2 has 
the value 4 J. In the 3D diagram, figure [T], this produces a kink at ^aver ~ 3.4 J. The 
phase boundary continues to grow, showing another slight change of shape around the 
minority component Van Hove singularity. In 2D (figure [2]), the critical value of 5/i has 
a pronounced maximum at the singularity at fiaver ~ 2.9 J, after which it decreases. The 
singularity in 2D is stronger than in 3D and the minority component does not reach 
its Van Hove singularity before half filling. In ID this also produces a clear feature at 
fJ'aver ~ 2 J, wherc the phase boundary essentially becomes horizontal, see figure [3l 

In our calculations we have taken the lattice height to be 2.5Eji and in the 2D 
and ID calculations we have taken the lattice height in the orthogonal directions to be 
10-E/j. We have used 6 for the mass number of the atoms, corresponding to Lithium. 
Our calculations are in the intermediate coupling regime, with U/J<6. We have also 
checked our results in the weak coupling BCS limit and found no qualitative difference 
to the results presented here. The exphcit coupling strengths {U / J) used are —3.7 in 
3D, —3.3 in 2D and —3.2 in ID. Note that all calculations are performed for the full 
3D system, just with the higher lattice heights in the orthogonal directions. Therefore 
the calculations take into account the small but finite tunneling between the one or two 
dimensional systems. 

3.2. Finite temperature 

We have studied the phase diagrams also in finite temperature, and our results indicate 
that the FFLO area, and features related to it, of the phase diagram, will gradually 
disappear with increasing temperature, as shown in figures HHSl In high temperatures 
even the effect of dimensionality seems to disappear, and phase diagrams b-d in[5]have no 
qualitative differences. It should be noted that in lattices the BCS critical temperature 
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Figure 4. Finite temperature phase diagrams in 3D and 2D lattices: (a) 3D at 
UbT/J w 0.1 (or 10 nK), (b) 3D at hgT / J « 0.2 (or 20 nK), (c) 2D at ksT/J « 0.1 (or 

10 nK), and (d) 3D at ksT/J « 0.2 (or 2()nK). The dashed Hne ( ) always shows 

the majority Van Hove singularity, where /i-f = 4J, and the dash-dotted line ( — • — ) 
shows the minority singularity, with /ij = 4,7 (only reached in 3D). The calculated 
data points at the phase boundaries are solid: BCS-FFLO and circles: FFLO - normal 
state. The temperatures in nK correspond to Jx = 0.07 Eji. 

Tc depends on the average filling factor or average chemical potential, since Tc is directly 
proportional to the gap at zero temperature. At half filling Tc is much higher than at 
low fillings. At finite temperatures, the BCS state can also support some polarization, 
due to thermal fiuctuations. 

4. Effect of the hcirmonic confinement 

We have used the local density approximation to study effects of a harmonic trapping 
potential superimposed on the lattice. Assuming that the harmonic trap is sufficiently 
shallow, it is possible to account for its effects by taking it to be locally constant 
and letting the chemical potential vary as a function of position, the standard local 
density approximation (LDA). This means that we take the system to consist of cTibical 
regions of N^^-^ lattice sites, where we approximate the local chemical potential for each 
component to be /u'^(r) = — V{r), where r is the distance from the center of the 
trap to the center of the cube, V{r) is the trapping potential, and ficr is the global 
chemical potential for component a, which is constant throughout the system. Because 
the system is approximated with an infinite lattice inside each region, the value of Ngrid 



FFLO state in 1, 2, and 3 dimensional optical lattices combined with a non-uniform background potentials 




Figure 5. Finite temperature phase diagrams: (a) ID at ksT/J « 0.2 (or 20 nK), 
(b) 3D at keT/J « 0.4 (or 40 nK), (c) 2D at keT / J « 0.4 (or 40 nK), and (d) ID 
at ksT/ J fa 0.4 (or 40 nK). In (a), the flat top (see figure [3]) disappears, because the 
maximum is so smaU that /i| does not reach the Van Hove singularity. In the high- 
temperature diagrams, (b)-(d), FFLO and all the features related to the dimensionality, 
have vanished. The temperatures in nK correspond to Jx = 0.07 En. 

can be chosen freely and does not have any physical implications. 

In this scheme, reading the effect of the trapping potential from the phase diagrams 
computed with respect fiaver and 6fi is exceedingly simple: because both of the 
components see the same trapping potential, 6fi does not depend on r, and any situation 
with given numbers of atoms in each component corresponds to a horizontal line in the 
appropriate phase diagram. Since the maximum number of identical fermions in the 
ground state of each lattice site is 1 for each component in the gas, the order parameter 
A goes to zero as the filling factor of either of the components approaches zero or one. 
This implies that with sufficiently high total particle numbers the region in the center 
of the trap is not superfluid, but a normal gas, surrounded by a superfluid shell, which 
is again surrounded by a shell of normal gas [33l|3lj. Because at zero temperature, with 
any filling besides close to an empty or a full lattice, FFLO is always between the BCS 
and the normal states in the phase diagram, there is a shell of FFLO between the BCS 
and normal regions in the trap. This is demonstrated in the gap profiles in figure El 
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Figure 6. Gap (A/J), average density ((ni +n2)/2), and density difference (ni —712) 
profiles for a 3D lattice in a harmonic trap with LDA, together with the FFLO wave 
vector q. The frequency of the harmonic trap is chosen to be 120 Hz and the wavelength 
of the lattice d, 515 nm. Shown arc four different shell structures: (a) with Sfi/J « 1.7 
and the density in the center of the trap at half filling, the system forms a polarized core 
in the FFLO state, surrounded by a shell of polarized normal gas. In (b), ~ 1.6, 

but there are more atoms overall in the system, giving a higher density in the center of 
the trap, which leads to a core of normal gas to form inside a shell of FFLO. Figure (c) 
shows an even more intriguing shell structure, corresponding to Sfj,/ J ~ 1.4: normal 
gas - FFLO - BCS - FFLO - normal gas. Note that this structure is dramatically 
reflected in the density difference. In (d), the system is unpolarizcd, but the density 
in the center of the trap is so high that a core of normal gas is formed inside a shell 
of BCS superfluid. All the plots correspond to horizontal lines in figure [TJ The overall 
polarizations P ^ {N^ - Ni)/{N^ + Ni) are: (a) 0.44, (b) 0.27, (c) 0.13, and (d) 0. 



5. Density-density correlations 

Here we compute the density-density correlation functions after free expansion. 
Correlation functions are defined by [2Zj 

Gn(r,r',t) = (nt(r,t)ni(r',t)) - {n^{r,t)){h^{r' ,t)) (5) 
= (^|(r,t)^^(r,t)^|(r',t)^i(r',t)) - {n^{r,t)){n^{r' ,t)), 

where f and | are the component indices and t is time. In the above formula we have 
subtracted the term with mean atomic densities, since often it is easier to focus on 
fluctuations. Because the density- density correlations after the free expansion reflect 
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correlations in momentum space at t = 0, correlations can be computed using the wave- 
function in momentum space prior to expansion [SH]. Positions in real space after time 
t are related to k- vectors trough r = htk/m. 

In momentum space it turns out that the plane-wave FFLO density-density 
correlations are strongly correlated at points which correspond to momenta k and 
— k -|- 2q, where q is the wavevector associated with the FFLO state. In the BCS 
state density- density correlations peak similarly, but with q = 0. The two-mode FFLO- 
state i.e. the state where the pairing gap is given by A(r) = Aocos(q ■ r), can also 
leave a clear signature on the density- density correlations. In the two-mode FFLO-state 
an f-atom which is in the momentum state k is paired with |-atoms which are in the 
momentum states — k + q and — k — q. This gives rise to large correlation peaks when 
k + k' ± q = 0. 

In figure [7] we demonstrate the difference between the BCS density-density 
correlations and the FFLO density-density correlations at zero temperature. In figure 
(a) we have plotted an example of the BCS-state density- density correlation between 
the components at the z = plane while in figure (b) we show on example of the 
FFLO density- density correlation between components, again at the same plane. As 
one can clearly see, the correlations are very different and the difference arises from 
two reasons. First, in the FFLO-state the density-density correlation peaks have been 
effectively shifted by 2q. Second, in the FFLO-state density-density correlation there 
are areas in the momentum space where correlation peaks vanish. The reason for this 
"breach" is that one quasiparticle dispersion has changed sign in the peakless region. 
Physically this means that these areas are populated only by normal atoms, and there 
are no pairs to give rise to correlation peaks. The height of the correlation peaks contains 
information on the underlying pair wave-function ~ ttk^k iSSj . In the weakly interacting 
BCS limit this function is strongly peaked at the Fermi surface, but it becomes broad 
and featureless in the BEC limit. In (c) and (d), we show that the difference between 
the BCS- and FFLO-state density- density correlations can persist even for integrated 
correlation signals 

C^,i{x,y) = J dz Gi^^i {x,y, z, —X + ht2qx/m, —y + ht2qy/m, —z + ht2qz/m) , (6) 

although integration obviously smooths out some features, especially the complete 
"breach" is not visible. At non-zero temperatures, sharp areas without correlation 
peaks disappear with increasing temperatures. However, the shift in the positions of 
the peaks persists even at non-zero temperatures |35j . 

The difference between density- density correlations of the FFLO-state and the BCS- 
state can be seen even more clearly when the lattice is one dimensional (ID) i.e. when 
tunneling strengths in y- and z-directions vanish. In figure |H] we show an example of 
density-density correlations in a ID lattice. There (a) demonstrates the antibunching in 
the BCS-state density- density correlation of a single component. The result is similar to 
the one in the ideal Fermi gas However, in the superfluid BCS state a "bunching" 
peak appears at Q{x) — Q{x') = 0, where Q{x) = mx/{ht) = k, where k is a lattice 
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momentum. This peak is absent for the ideal Fermi gas. In (b), we demonstrate 
the BCS-state density-density correlation between the components and in (c) we show 
the FFLO-state density- density correlation between the components. As one can see 
from figure (b) the the BCS-state density-density correlation between components is 
symmetric with respect to x = and it has no regions without correlation peaks. On the 
other hand, from figure (c) we can see that the FFLO-state density-density correlation 
has a region (marked with dashed lines) without correlation peaks and one can also see 
that the FFLO-state density-density correlation is not symmetric with respect to x = 0. 
Note especially that in ID the region empty from correlation peaks is not vanished by 
integration over other dimensions. In figure (d) we show that this remarkable signature 
appears in the low filling well, which corresponds to the continuum ID system 

that is of current experimental interest as well. 

Interestingly, the phase-separation between the normal gas and a paired-state could 
be visible in the density-density correlations between components. This follows from 
the fact that the density-density correlations between components in the normal state 
vanish, whereas in the paired state the correlations are at their strongest around the 
Fermi momentum. In a lattice superimposed by a trap the local density, and therefore 
the local Fermi momentum, is different in different areas of the gas. This may allow to 
identify spatial phase separation and shell structures of normal and paired states from 
the freely expanded cloud where momentum has been mapped into position. 



6. Fixed density phase diagrams 



In ultracold gas experiments, the number of particles is often the fixed quantity, not the 
chemical potential. Although the trapping potential makes also fixed chemical potential 
calculation relevant in the LDA sense (as discussed in section Hj), it is of interest to 
consider also the fixed particle number case, because in the case of optical lattices the 
background potential is a practical issue which may be eliminated to a certain extent. We 
have studied the phase diagrams as a function of polarization P = (n^ — ni)/{n-i + 
and temperature, i.e. keeping the filling factors rij and constant with respect to 
position. 

In this situation, in addition to minimizing the free energy F, it is necessary to 
solve the chemical potentials from the number equations: 
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Figure 7. Differences between the BCS-state density-density correlations and the 
FFLO-state density-density correlations in a 3D lattice. In (a) is shown an example 
of the BCS-state density-density correlations in the z = plane, while in (b) we show 
an example of the FFLO-state density-density correlation in 2 = plane. In (c) and 
(d) are shown the corresponding signals which have been integrated over z. In (a) 
and (c) we chose r — — r' and in (b) and (c) r = — r' + Ht2q/m. Note that any 
other choice would produce zero correlation. All these examples were calculated at 
zero temperature and in (a) and (c) P = 0.0, (n| -|- ni)/2 = 0.55, and A/(2J) — 0.49 
and in (b) and (d) P = 0.168, {n^ + n|)/2 = 0.55, 2q^ = 0.25(7r/d), qy = = 0, 
and A/(2J) = 0.16. Color-coding is such that warm colors imply high peaks and cold 
colors low, but in the white areas the correlations vanish identically. 

and / is the Fermi function. This scheme produces the following stable phases: a 
polarized superfluid phase, which is the standard BCS phase when P = and the 
BP/Sarma phase when P > 0, FFLO, and normal state. It is important to note that 
when the densities are kept fixed, the BCS state is a special case of the BP state, with 
zero polarization. In addition to these phases, we have considered a phase separated 
state consisting of an unpolarized BCS gas and a polarized normal gas, as suggested 
in [36]. Such phase diagrams have already been discussed in [23] for a 3D lattice, here 
we expand those results to a ID system and discuss the effects of Hartree and Gorkov 
corrections on the phase diagrams. 

The temperature-polarization phase diagram in a ID lattice is shown in figure 
ini The interaction is chosen so that the BCS critical temperature is close to that in 
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Figure 8. ID lattice density-density correlations. In (a) we demonstrate the 
antibunching effects of the BCS density-density correlation of a single component 
in a ID lattice. In (b) is shown the BCS-state density-density correlation 
between components and in (c), we show the FFLO-state density-density correlation 
between components and (d) the FFLO-state density-density correlation between the 
components with low average filling fraction. In (c) and (d) dotted lines indicate 
the gapless regions. In (a) and (b) the polarization P = 0, (ri| + ni)/2 = 0.41 and 
A/(2J) = 0.53. In (c) the polarization P = 0.48, (n| + ni)/2 = 0.40, 2q^ = 0.40(7r/d) 
and A/(2J) = 0.19. In (d) the polarization P = 0.91, (n| + ni)/2 = 0.20, 
2qx = 0.37(7r/(i) and A/(2J) — 0.041. In (a), Q{x) = mx/{ht) and we have chosen 
X = 0. In (b), we have chosen x + x' = and in (c) and (d) we have chosen 
X + x' — 2htqx/'rn = 0. For other choices, correlations vanish. The distance between 
the peaks in figures 8b-d, as well as in figure 7, reflect the discreteness of the finite size 
lattice. 
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the diagrams for the 3D lattice, shown in [TOl These figures demonstrate that with 
comparable conditions, a ID lattice can support much higher critical polarizations in 
FFLO and polarized superfiuid phases than the 3D system (~ 0.8 vs. ~ 0.3 in zero 
temperature with the chosen parameters). 

The Hartree corrections rise from the Hartree terms, f^(c|i^C|k)c|i^C|k and 
^c|,^C|k(c|,^C|k), which have been left out of the Hamiltonian. When the density is 
constant in the system, the effect of the Hartree terms is implicitly included in the 
number equations, (171). However, as the phase separated state contains components with 
different densities |23], including the Hartree terms explicitly may change the difference 
in free energies between phase separation and e.g. FFLO. The free energy of system 
with constant densities and n^, with the Hartree terms included, is 

F = - — + /i^ni + fi^ni + 2Un^ni + 

E {^i-^^- + ^-••^.'^ - ^ + ''''""'"^ + e/^^— ))) (8) 

We have found that while including the Hartree corrections brings the absolute 
difference, Fps — Ffflo closer to zero, it enlarges the FFLO area in the phase diagram, 
compared to phase separation. This is shown in figure [TUl where two phase diagrams 
with identical conditions, except for the inclusion of the Hartree terms, for a 3D lattice, 
are shown. 

For a two-component Fermi gas with equal densities, the transition temperature 
Tco to superfiuid state has been calculated as 

T,oOcEpexp(l/iV(0)f/o), (9) 

where iV(0) is density of states at Fermi surface and Uq is the two-body interaction. 
This transition temperature is obtained by considering a purely two-body interaction. 
However this picture can be modified due to the effect of the medium. This effect, 
which is also referred to as induced interaction [37j, was originally studied in 60s by 
Gorkov and Melik-Barkhudarov [3H]- The main idea goes back to the polarization in 
the medium by one fermion and its influence on a another atom. The second atom is 
scattered by a modifled total interaction, including induced interaction, which can be 
written as 

Utot = Uo + Uind. (10) 

Gorkov et al. found that the contribution of this effect reduces the transition 
temperature by a factor 2.2. In our formalism, the interaction enters via the flrst term 
in the grand potential |H Therefore in principal one could take into account the induced 
interaction correction by considering Utot instead of Uq. This requires to calculate the 
induced interaction which is ~ Uq and involves the Lindhard function. 
The induced interaction is [3D] 

^.n. = Ulm = Ul [ ,,^^,3^'~^'r (1^) 

J [2TinY ep+k - Cp 
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Let us here qualitatively discuss the effect having a lattice dispersion ek instead of the 
usual quadratic dispersion in free space. In the low density limit, the lattice dispersion 
becomes effectively quadratic, and one can expect the same factor of 2.2 reduction 
of critical temperature as in the free space case [3H]- The term in the denominator 
of the Lindhard function is, in case of the usual quadratic dispersion, of the form 
(k + p)^/2 — p^/2 = k'^/2 + kpcos6 (where 6 is the angle between k and p). Around the 
point ki = 7t/2 (for ID lattice at half filling), the lattice dispersions 1 — cos{ki) (here i 
means x,y,z depending on the dimensionality) become effectively linear. Linearizing the 
lattice dispersion produces to the denominator of the Lindhard function terms of the 
form \ ki+pi\ — \pi\ = a/ kf + pf + 2kiPi— \pi\. Assuming that this can be approximated by 
a Taylor expansion, one ends up with k"^ /2-\-kpcos6 like in the quadratic dispersion case. 
This indicates that the linear dispersion regimes might not change the Gorkov correction 
considerably. However, only a numerical evaluation of the integrals can give a definite 
answer. An interesting issue is the behaviour of the L{k) when the lattice filling is high 
enough to include in L{k) dispersions near the band edges, i.e. VanHove singularities. 
Then, the dispersion 1 — cos(fcj) around ki = n becomes quadratic again, but with 
a negative effective mass. Naively, this means the induced interactions would have 
opposite sign than usual, i.e. enhance Tc rather than suppress. However, the Lindhard 
function in this case, in addition to quadratic dispersion with negative mass, contains 
also additional terms due to expansion around vr, and again a numerical approach should 
be applied. Anyhow, one could still expect the Van Hove singularity to reduce the effect 
of the Gorkov correction, as has been shown in |10] for a 3D gas in a ID lattice. Note 
that this is not the same system as what we mean here by ID lattice, (which is a ID 
gas in a ID lattice) as explained in the introduction. In our ID lattice, the Van Hove 
singularity would not influence the Gorkov correction since the singularity happens at 
full filling, and that should be the same as low density limit due to symmetry of pairing 
of particles and holes over half filling. However, in 2 and 3 dimensions the Van Hove 
singularities as well as regions of linear dispersion affect L{k) and it is an intriguing 
question how these effects add up. Numerical investigation of these issues is a topic of 
our further work. 

7. Conclusions 

We have considered phase diagrams of density imbalanced two-component Fermi gases in 
optical lattices of dimension 1, 2 and 3. The phase diagrams in the plane of the average 
chemical potential of the two components, and the chemical potential difference, show 
striking effects originating from the VanHove singularities of the lattice. These features 
appear only for the FFLO state, not the BCS. Therefore they are unique signatures of 
the FFLO state and reflect the fact that the nesting of the Fermi surfaces in a lattice 
enhances FFLO pairing compared to case of the homogeneous space. We show how these 
features preserve to finite temperatures and finally disappear for very high temperatures. 
Using LDA, we have demonstrated various shell structures that can appear when the 
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T(nK) 

Figure 9. The phase diagram of an imbalanced Fermi gas in an effectively ID lattice. 
Colors: polarized superfluid = blue (this is BCS at P = 0), FFLO = yellow, phase 
separation (PS) = red, normal — white. The average filling factor is 0.2 atoms/lattice 
site in each component, Jx = Q.OTEn, and U = — 0.2£'fl;, where Eji — H^k"^ /2m is the 
recoil energy. The FFLO area is remarkably large, with a high critical polarization. 




10 20 30 40 10 20 30 40 

T (nK) T (nK) 



Figure 10. Phase diagrams in a 3D lattice, without (left) and with (right) the Hartree 
corrections included. The colors are the same as in figure [5J polarized superfluid = 
blue, FFLO = yellow, phase separation (PS) = red, normal = white. The average 
filling is 0.2 atoms/lattice site in each component, J — Q.QlEn, and U = — 0.26£'fl;. 
The phase diagram on the left is published in f23' and is shown here for comparison 
to the one on the right. 

lattice is superimposed by a harmonic trapping potential. For the studies of the FFLO 
state itself, it is useful that structures where only FFLO and normal states appear can be 
found, but also more exotic shell structures such as normal-FFLO-BCS-FFLO-normal 
are possible. 

Density-density correlations are one possibility for observing the various phases and 
states. We have shown here how the unpaired atoms in the FFLO state leave a clear 
signature in the correlations. Especially in ID the signature is very prominent. This is 
true for the ID gas in a ID lattice as considered here, as well as for the ID continuum 
system (our system in low density limit) which is also of high interest since the one- 
dimensional confinement is known to enhance FFLO pairing [HI |l2] even without the 
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lattice potential. 

We also considered Hartree corrections which according to our calculations tend 
to increase the FFLO area versus phase separation in the case of fixed particle 
numbers. Outlook for further work inchides more detailed considerations of the 
Gorkov corrections, and studies of strongly interacting gases instead of the weak and 
intermediate coupling regime considered here. 
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